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Abstract. We consider the six- vertex model on an N x N square lattice with the 
domain wall boundary conditions. Boundary one-point correlation functions of the 
model are expressed as determinants of N x N matrices, generalizing the known result 
for the partition function. In the free fermion case the explicit answers are obtained. 
The introduced correlation functions are closely related to the problem of enumeration 
of alternating sign matrices and domino tilings. 
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1. Introduction 



The six-vertex model was studied for both periodic [I], |2| and fixed boundary conditions 
@- B B 0- ^he particular example of fixed boundary conditions J7j of the model on 
aniV x JV square lattice is the so-called domain wall boundary conditions (DWBC) M. 
Under special restrictions on the vertex weights this model is related to the enumeration 
of alternating sign matrices [H| and domino tilings of Aztec diamonds |§. 

The model with DWBC originally appeared in the context of investigation of 
norms of the Bethe states in the framework of the Quantum Inverse Scattering Method 
(QISM) pTT[| . In the last decade the six- vertex model with DWBC has found interesting 
applications in different fields of physics and mathematics [fT2 , |13| , |H] , [15|, [T(| due to 



the results of the papers [17], [T8[], where the determinant formula for the partition 
function has been obtained and proved. This determinant formula allowed to solve 



several problems in combinatorics pi which were standing for a long time [19 



A wide range of problems such as refined enumeration of alternating sign matrices 
(ASM) 0] and the arctic circle theorem |21], ^] can be solved only if the correlation 



functions of the model are known. In general, the calculation of the correlation functions 
is a more complicated problem than that of the partition function. Additional difficulties 
may arise due to the lack of translation invariance caused by the fixed boundary 
conditions. 

In this paper we will consider two kinds of one-point boundary correlation functions 



of the six-vertex model with DWBC. The function of the first kind, G 



(M) 
N 



is the 



local state probability on the boundary vertical edge and it may be called "boundary 
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spontaneous polarization". The function of the second kind, , describes the 

probability of the vertex being in the specific state. For the model on an N x N 
square lattice we obtain representations for these correlation functions as determinants 
oi N x N matrices. These correlation functions are the generalization of a boundary 



correlation function considered in [23 



There are three convenient ways for description of the six-vertex model: (i) in terms 
of arrows pointing into and away from each vertex; (ii) in terms of lines flowing through 
the vertices; (iii) in terms of spins on the edges. The six types of vertices allowed in the 
model are plotted in Figure [I]. A statistical (vertex) weight corresponds to each type of 
vertex. We consider the six- vertex model with the vertex weights being invariant under 
the simultaneous reversal of all arrows. Hence, there are three different vertex weights, 
a, b and c, see Figure 1. 

The domain wall boundary conditions imply that all arrows on the top and bottom 
of the lattice are pointing inward while all arrows on the left and right boundaries are 
pointing outward, see Figure @(a). It means that the solid lines flow from the top of 
the lattice to the left boundary, see Figure |](b). In Figures ||(a) and ||(b) one of the 
possible arragements is presented as well. The domain wall boundary conditions in 
terms of spins are shown in Figure 0(c). 

In the inhomogeneous model the wertex weights a, b and c are site dependent. To 
introduce this dependence we will use two sets of the variables {X a } and {^} that are in 
one to one correspondence with the set of lines. The rows will be enumerated by Greek 
indices a = 1, . . . , N and the variable X a corresponds to a-th row; the columns will 
be enumerated by Latin indices k = 1, . . . , N and the variable corresponds to fc-th 
column. This correspondence is shown in Figure 0(b). Each statistical weight associated 
with the vertex lying at the intersection of a-th row and /c-th column will depend on 
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Figure 1. The six allowed types of vertices: in terms of arrows (first row), in terms 
of lines (second row) and in terms of spins (third row) . 
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Figure 2. One of the possible configurations in the model with DWBC: 
of arrows; (b) in terms of lines, (c) DWBC in terms of spins. 



(a) in terms 



the pair of variables (A Q , u k ). The parametrization that allows one to apply QISM is 
a(A a , u k ) = sinh(A Q - u k + rj), 

b(A Q , u k ) = sinh(A Q - u k - rj), (1) 
c(Aq,, u k ) = sinh 2rj. 

In the homogeneous limit all X a — > A, and all u k — > u. All positive values of the vertex 
weights, up to an overall scaling transformation, may be obtained by choosing both 
A — u and r\ either real or pure imaginary. 

Our calculations are based on the Quantum Inverse Scattering Method (QISM), 
reviewed briefly in Section 2. In Section 3 we derive the reduction formulae for 
the boundary correlation functions G^P and The recursion relation for the 

partition function follows from these formulae as a particular case. The determinant 
representation for G^ and is obtained from these reduction formulae in Section 4. 
In the free fermion case the homogeneous limit for these boundary correlation functions 
is calculated explicitly in Section 5, while the general case of the homogeneous limit is 
considered in Section 6. 



2. Formulation of the model within QISM formalism 

To apply the Quantum Inverse Scattering Method [[l l| we use the spin description of 
the model. With each vertical line (column) and horizontal line (row) one associates 
the space C 2 , with spin up and spin down states forming a natural basis in this space. 
The total space of the vertical lines is V = (C 2 )® N and the total space of the horizontal 
lines is H = (C 2 )® N . With each vertex of the lattice one associates an operator acting 
in the full space V <8> Ji. This operator is called L-operator and it acts nontrivially only 
in a single horizontal space C 2 and in a single vertical space C 2 , while in all other spaces 
it acts as the unity operator. To distinguish the spaces in which the L-operator acts 
nontrivially one can label it as L ak and associate it with the vertex being the intersection 
of a>th row and k-th column. The matrix elements of the L-operator (which is 2 2N x 2 2JV 
matrix) are either zeros or functions a(A Q , u k ), b(A a , u k ), c(A a , u k ), defined in ([]]). Hence, 
the L-operator L ak is the function of \ a and u k , L a j fc (A a , u k ). Since the L-operator acts 
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nontrivially only in the direct product of a pair of two-dimensional spaces, all its elements 
may be written in the compact form as 


b(\ a ,v k ) c(X a ,u k ) 
c (A a , v k ) b (A a , i/ fc ) 




L Q fc(A Q , Uk) 



/a (A Q , v k ) 












(2) 



(X a ,v k )/ [afc] 



This is the matrix with respect to a-th copy of C 2 in 7i and fe-th copy in V, with the 
matrix elements being trivial matrices in the rest copies of C 2 in V and in Ji. One 
can write the L-operator in the alternative form with the separated "horizontal" and 
"vertical" spaces, namely, as a matrix with respect to a-th copy of C 2 in 7i with the 
operator matrix elements acting nontrivially only in fc-th copy of C 2 in V 

sinh(A Q — v *. + r\o^) a k sinh 2r] 

a£ sinh 2rj sinh(Ao, — v k — rial. 



L Q fc(A Q ,, Vk) 



(3) 



where cr|, a k — \{p% ± ierf) are Pauli matrices. 

The main object of QISM is the "vertical" monodromy matrix T a (A a ) which is 
defined as the ordered matrix product of the L-operators along a-th horizontal line 



T Q ,(A Q! ) — Lo,jv(A a , vn) ■ ■ ■ L Ql (A a , vi) — ^ j j ^ 



(4) 



All entries of the monodromy matrix T Q ,(A Q ) are operators acting in V and they depend 
on the variables v\, . . . , vn, i.e., A(X) = A(X; {vk}k=i)i e ^ c - Obviously, instead of the 
"vertical" monodromy matrix one may use the "horizontal" one, which is the ordered 
product of the L-operators along the vertical line. 

The Quantum Inverse Scattering Method is based on the intertwining relation for 
the L-operators: 

Rap(X a , Xp) L Q fc(A Q , Vk) L/3k(X/3, v^ = L/3k(X/3, Vk) L Q fc(A Q , Vk) R a p(X a , Xp), a ^ (3. (5) 

The R-matrix R a p(X,X') acts nontrivially in the direct product of a-th and /3-th 
horizontal spaces and are given by 

//(A', A) 



Rq/9(A, A') 



V 










1 

<?(A',A) 




where the functions /(A', A) and g(X', X) are 
/( A,A) = -"(A-A< + 2„) 





<?(A',A) 
1 




i?(A',A) 









\ 



(6) 



/(A', A) J M 

sinh 2rj 
' sinh(A - A') ' 



sinh(A-A') ' ^"''v sin h(A-A')' ^ 

This is the so-called trigonometric R-matrix [11|, which satisfies the Yang-Baxter 
equation 

Rap(X a , Xp) Ro 7 (A Q , A 7 ) Rpj(Xp, A 7 ) 

= Rp~f{Xp, A 7 ) R Q7 (A Q , A 7 ) R a p(X a , Xp), a^^7- (8) 
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Due to relation ([5]) and commutativity of the matrix elements of L-operator (|3]) at 
different lattice sites one has the intertwining relation for the monodromy matrix: 

Rap(^a, A/?) T Q (A Q ) T^A^) = T^A^) J a (X a ) Rq-/3(A q , A/?), a ^ (3. (9) 

Equation (0) defines the commutation relations for the operators entering the 



monodromy matrix. The complete list of these relations can be found, e.g., in (TT 
For our purposes we need only two of them: 

A(X) B(X') = /(A, A') B(X') A(X) + g(X', A) B(X) A(X'), 

B{X)B{X') = B{X')B{X). 

As generating vector in the space V it is convenient to use the state either with all spins 
up or with all spins down 

i K) = <=ii n* = ®ti (J) , m = <=ii i )* = ®£i (°) ^ • (ii) 

These vectors are annihilated by the operators C(A) and B(X), respectively, 

C(A)|t> = 0, S(A)|^) = 0, (12) 
and they are eigenvectors of the operators A(X) and D(X) 

A(X)\f) = a(X)\f), D(X)\t} = d(X)\t), 
A(\)\^) = d(\)\D,), D(A)|^) = a(A)U), 
where the functions a(A) and d(A) are equal to 

N N 

a(A) = JJsinh(A - i/ fe + 77), d(A) = sinh(A - v h - 77). (14) 

k=l k=l 



(13) 



The vectors ( ff | and ( Jj |, dual to ([TT]) , are eigenvectors of the operators A(X) and -D(A) 
with the same eigenvalues as in equations ( |T3"D while instead of equations (|l"2"D one has 

(ft | 5(A) =0, (^|C(A) = 0. (15) 

Consider vectors generated by multiple action of operators B(X a ) on the state Iff) 

73(Am)...5(A 1 )|^), M<iV. (16) 



The result of the action of the operator A(X) on vector (16 ) follows from commutation 
relations fllOD 



M M M M 

a(x) n B(x a ) 1 1 > = a n B(x a ) a ^( a ) n b ^ \ ^ w 

a=l a=l 3=1 a=1 

where the coefficients A and hp are 

M M 

A = a(A) J] /(A, A 7 ), A^ = a(A /3 ) g{Xp, A) J] /(A„, A 7 ). (18) 

7=1 7=1 

The partition function Z N = Z N (Xi, . . . , A^; ui, . . . , i^v) of the model on an N x N 
square lattice with DWBC is obtained by summation over the contributions of all 
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possible spin configurations. The contribution of each configuration is equal to the 
product of all vertex weights of this configuration. In terms of QISM the partition 
function may be represented as 



§ «<T| J ® ( S k(i I J T N (X N ) . . . Ti(Ai) ( ® | T>* J ® ( ® U) a ) • (19) 

a=l / \k=l I \k=l I \a=l 



The boundary conditions on the left and the right boundaries extract from each matrix 
T a (A Q ) the operator B(X a ). The boundary conditions on the top (bottom) of the lattice 
correspond to the vector | ff ) (( JJ. |), respectively. Hence, the partition function can be 
written in the form 

Z N = (^\B(X N )...B(X 1 )\f). (20) 

Due to relation flip]) the order of operators B(X a ) in the product is not essential. 

The determinant representation for the partition function was obtained in the 
papers [T7|, |H| and has the form 

N N 

f] EI sinh(A Q -v k + v) sinh(A Q - u k - rj) 

Z n = a=1 JL =1 — — — detTv^. (21) 

[[ smh(A /3 - A a ) 11 smh(v k - uj) 

l<a</3<N l<k<j<N 

The entries of the matrix Z are given by 

Z ak = (j)(X a ,v k ), a,k — l,...,N (22) 

where the function 0(A, v) is 

x sinh2r? 

^ v ) = — u\ Z S • a v ( 23 ) 

smh(A — v + rj) smn(A — v — rj) 

The proof of determinant representation ( pi]) based exclusively on commutation relations 
(p!0|) is given in Section 4. 



3. Boundary correlation functions 

In the present paper we consider two kinds of correlation functions describing the 
local state probabilities at the boundary. The first correlation function describes the 
probability of absence of vertical solid line between M + 1-th and M-th rows on the first 
column and is known as "boundary spontaneous polarization". In terms of QISM it is 
the one-point correlation function of the local spin projector q 1 = |(1 — erf) on the spin 
down state, and it can be written as 

G ( P = Z N l (ij-\ B(\ N ) . . . B{A M+1 ) Ql B(\ M ) . . . B(X 1 ) (24) 

The second correlation function describes the probability that the solid line on the first 
column turns to the left just on M-th row 

H^ I] = Z N 1 (\^\ B(X N ) . . . B(X M+ i) qi B(X M ) Pi £(A m _ x ) . . . S(Ai) (25) 
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where p\ is the projector on the spin up state, p\ = |(1 + erf). Since pi + q± = I, these 
correlation functions are related to each other as follows 

G! v M) = < ) + <- 1) + ... + ^ 1) , (26) 

rr(M) _ r (M) r (M-l) (e)7 s 



However, it is easier to calculate them from definitions (|24| ) and ( p5|) independently. 

In this Section we express the correlation functions of the model on an N x N square 
lattice through the sum over partition functions of the models on (N — 1) x (N — 1) 
square sublattices. We will call the corresponding formulae the "reduction formulae". 
The derivation of these formulae is based exclusively on commutation relations fllOl) . 
Since G^ = 1, in the particular case M = N the reduction formula for turns 
into the recursion relation for the partition function. In the next Section we will 
prove that determinant representation ( pT|) is the solution of this recursion relation. 



The determinant representation for the correlation functions can be obtained then 



by substituting expression (]2l|) in the reduction formulae, what makes our algebraic 
approach self-containing. 

To derive the reduction formulae for the correlation functions we rewrite them 
in the form suitable for applying commutation relations ([K]). Let us decompose the 
monodromy matrix T a (X a ) into the matrix product of two monodromy matrices in a-th 
space 

T a (X a ) = T a2 (A a )T Q , 1 (A Q ,). (28) 

This decomposition of the monodromy matrix is known as the "two-site model" |^TJ. In 
our case we choose these matrices defined as follows 

T Q 2(A a ) = \- aN (X a , v N ) . . . L a2 (A 

VC 2 (A a J V 2 {A a )J [a] 
J (X)-l (X u)-( MK) Bl(K) ) 

The matrix elements of T Q 2(A) commute with the matrix elements of T a i(A) since 
they are operators that act nontrivially in different spaces. The entries of T Q , 1 (A) act 
nontrivially in the first "vertical" space and depend on the "vertical" variable Vx, while 
the entries of T a2 (A) act nontrivially in the rest N — 1 "vertical" spaces and depend 
on the "vertical" variables z/ 2 , . . . , z/jy. Each set of operators entering the monodromy 
matrices satisfy the commutation relations given by (f|) and, in particular, (10). It 



should be noted that the generating vectors | ft) ( | JJ- ) ) can be represented as the direct 
product of two generating vectors, e.g., | it) = | IT2 ) ® | fTi ), where | f|~2 ) = ®fcL 2 | T )k 
and | fix ) = | | )i- The properties of the both states | ff 2 ) ( | JJ-2 )) an d | 1Ti ) ( I JJ-i )) are 
the same as of | ft) ( | JJ-)), see the previous Section. The eigenvalue of, e.g., operator 
A 2 (X) on the state | -fo ) is 

TV 

A 2 (A)|ft 2 ) = a 2 (A)|ft 2 ), a 2 (A) = J] sinh(A - z/ fc + ??). (30) 

fc=2 
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The decomposion of the operators -4(A), -B(A), C(A) and -D(A) follows from definitions 
(p8|), (p^) and (f|). In particular, one has 

B(A) = A 2 (X)B 1 (X) + B 2 (X)D 1 (X). (31) 

The operators -Bi(A) and -Di(A) are the corresponding entries of L-operator (§) and they 
are given by 

£?i(A) = &l sinh 2% D x ( A) = sinh (A - z/i - 77^) . (32) 

Using formulae fl3~T|) and ( |32|) one can reduce the problem of calculation of the scalar 
products in the right hand sides of expressions (|24]) and ( |25|) to the problem of calculation 
of the scalar products involving the operators A 2 (X) and B 2 (X) only. Now we are ready 
to derive the reduction formulae for the correlation functions. 

We start with the derivation of the reduction formula for the function since 
this derivation is straightforward. Let us substitute expression ( |3"TD in (^) and calculate 
the scalar product with respect to ( 4j-i | and | ffi ) using formula (|32|) . We are left with 
the expression 

M-l N 

H^ 1 = Z^ 1 sinh2?7 J^J sinh(A Q — v\ — if jQ sinh(A Q — v\ + if 

a=l a =M+l 

x (h\B 2 (\ N )...B 2 (\ M+1 )A 2 (\ M )B 2 (\ M _ 1 )...B 2 (\ 1 )\ii 2 ). (33) 

Applying (|T7| ) and taking into account (|30|), we reduce the scalar product in (|33|) to 
the sum over scalar products that involve only the operators B 2 . Comparing these 
scalar products with expression fl20|), one immediately gets the following formula for the 
correlation function H N 



(AO. 
N ■ 

M-l N 



= Z^ 1 sinh 2i] J^J sinh(A Q — v\ — r/) J J sinh(A Q — ui + r/) 

a=l Q =M+1 
M <\ \ \ M 

x £ ! ?' x II A ->) ^-1 ({AJl W/3 ; 0*}£J • (34) 

/(A/3, Am) ^ 

In derivation of this formula we have used the fact that the functions /(A', A) and g(A', A) 
defined in (^) satisfy the condition g(X, A) //(A, A) = 1. This allows us to include the 
first term arising in (O) into the sum over (3. 



Equation ( 34] ) expresses the particular boundary correlation function of the model 
on an N x N square lattice as the sum over partition functions of the models on 
(N — 1) x (N — 1) square sublattices. It should be stressed that in derivation of 
the formula we have used only the algebra of operators 4 2 (A) and B 2 (X) given by 
commutation relations (|10|). 

Let us turn now to the correlation function G^ 1 ' defined by expression (^4j). At 
first, we substitute formulae fl3~T|) and (^) in (p4|). This gives 
M p-i N 
= Z^f sinh 2r] Y[ sinh(Ao, — v\ — vf j j sinh(A a — v\ + rf) 

/3=1 a=l a=/3+l 

x (^2 I B 2 (X N ) . . . B 2 (X P+1 ) A 2 (X P ) B 2 (X^x) . . . B 2 (\i) | ft 2 ). (35) 



Boundary correlation functions of the six-vertex model 9 

To obtain the representation of the correlation function G^ as the sum over the 
partition functions of the models on (N — 1) x (TV — 1) square sublattices one can 
substitute formula in expression fl35|), but then the double sum will appear in the 
resulting expression. To represent the right hand side of expression (|35|) as a single sum, 



like in formula (|34|), we choose another way. In expression fl35|) we pick up the term 
which contains the operator A 2 standing outmost to the left among the operators B 2 , 
i.e., the term corresponding to (3 = M. Taking into account formula ( |T7| ) it is easy to 
see that this term contains the following contribution, which corresponds to the first 
term in the right hand side of (|Hj): 

( ^2 I B 2 (X N ) . . . B 2 (X M+1 ) A 2 (X M ) 5 2 (A M -i) • • • B 2 (X 1 ) | ft 2 ) 

A/-1 

= a 2 (X M ) 11 /(A M , A 7 ) ( ^ 2 | B 2 (X N ) . . . B 2 (X M+1 ) B 2 (X M ^) . . . B 2 (X 1 ) \ ft 2 > 
7=1 

+other terms, (36) 

where the "other terms" are the terms with scalar products involving the operator 
B 2 (Xm)- Substituting equation ( |3~6"1) in ([35|) and applying formula (|17D to the rest terms 
in (|35|) , we obtain that 

M N 

= Z^ 1 Y[ sinh(A a - v x - 77) sinh(A a -1/1+77) 



a=l a=M+l 
M-l 

I a 2 [ A m ) sum 2r? 

X 



a 2 (A M ) smh 277 n 
smh(A M -Vi-rj) ±± 



x(h\ B 2 (X N ) . . . B 2 (X M+1 ) B 2 (X M -i) ■ ■ ■ B 2 (X 1 ) \f 2 ) + other terms ) , (37) 



where the "other terms" are again the terms, different to those in (|36|), with scalar 
products involving the operator B 2 (Xm)- Thus, the contribution written explicitly 
in (j3~7l ) is the only possible one which does not contain the operator B 2 (Xm) in the 
scalar product. Obviously, this contribution is symmetric under the permutations of 
the elements in the set Ai, . . . , Xm-i- At the same time the correlation function 
is symmetric under the permutations of the elements in the set A 1; . . . , A^, due to the 
commutativity of the operators B 2 (X a ). It follows immediately from these symmetry 
considerations that the whole expression standing in the parentheses in (|37| ) is just the 
sum over the cyclic permutations of the elements in the set Ai, . . . , Xm of the term that 
written down explicitly in the right hand side of equation (p7[). Therefore, the following 
representation for the correlation function G^ is valid 

M N 

G^ ] = Z N Y ] J sinh(A Q -v x — rj) \ j sinh(A a - v\ + 77) 

a=l a=M+l 

X sinUX^ - v -n) A f^ 13 ' X ^ Zn - 1 {{ X <*}a=l,a^ {Vk}k=2) ■ ( 38 ) 
3=1 ^ 13 1 7=1 

1*0 
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The procedure described above is similar to the derivation of equation fllTD from 
commutation relations fllOD , see, e.g., [|TT|. 



The important point to be discussed now is the meaning of representation (|3"8"D at 
M = N. Since the left hand side of equation (^) is equal to one at M = N, it turns 
into the recursion relation on the partition function. In this relation is expressed 
through the sum over the partition functions Z N _i 

N N N 

Z N ({\ a }a =1 ] {vk}k=i) = sinh2?7^ sinh(A a -v x - rj) sinh(A/3 -Vk+v) 



8=1 a=1 fc=2 



N 

X 



I I /(A/,, A 7 ) x Z N _, ({A a }? =W/3 ; K}f =2 ) ■ (39) 



7=1 
7^/3 



This formula was obtained by slightly different technique in |[24|| . Due to the symmetry 
of Zn with respect to the permutations of the variables {vk}k=i one can rewrite this 
recursion relation in a general form 

N N N 

z n ({A a }a=i5 { u k}k=i) = sinh 2rj ^ Yl sinh (^ ~ v j ~ v) sinh(A /3 -v k + v) 



8=1 a=1 fc=1 



N 

X 



J] /(A/3, A 7 ) x Z N -i ({AJl W/3 ; R}f=w;) , j = 1 , . . . , iV. (40) 



7=1 
7^/3 



Relations fl39|) and d4H|) are valid for arbitrary values of the variables {A a }^ =1 and 
{^fe}fcLi- Iterating relation fliPP iV — 1 times, with the initial condition Zi = sinh2?7, 
one gets the answer as the sum over all permutations P of {X a }a=i (c.f., (3.17) of J7[) 

Zjy = (sinh277) Jv y^ j [ f(Xp a , \p p ) sinh(A P(3 - u a - rj) sinh(A Pa - vp + rj). (41) 

P l<a<8<N 

On the other hand, one can show that determinant representation Q2T| ) for the partition 
function is the solution of recursion relation fl4T)|). The determinant solution of equation 
(HO) leads to the determinant representations for the boundary correlation functions. 

4. Determinant representations for boundary correlation functions 



We begin this Section with the proof that determinant formula (|2l|) for the partition 
function solves recursion relation (HO). The determinant representation for the boundary 



correlation functions will follow then from reduction formulae (p4|) and (38), which can 



be viewed as generalizations of the recursion relation for the partition function. 



Let us give the proof for j = 1 in (£5]). This case is given by relation (|39|). For 
the other values of j the proof is essentially the same. Consider the left hand side of 
relation (^) which is given by formula (|2T|). The aim is to represent this expression in 
the form given by the right hand side of relation fl3"9"|). Since Z N and Z^-i are expressed 
through determinants of N x N and (N — 1) x (N — 1) matrices respectively, it is clear 
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that relation (|39| ) can be proved developing a determinant of some N x N matrix by 
the first column. To apply this approach one has to transform the determinant of the 
matrix Z in ([21]) first. Consider the function 



N 

I] sinh(A 7 - A + 277) 
7=1 

N 

n sinh(A - v k - 77) 

k=l 



(42) 



For this function one has the identity 



N 



9n 



(A Q ) = J^$ fc 0(A c 



a 



,N, 



(43) 



k=l 



where the function 0(A, v) is given by formula (|23|). Here the coefficients are 
independent of a and are given as 



N 



sinh(A 7 - v k + 77) 

7=1 



A' 



(44) 



n sinh(z/ fc - z/j 



3=1 



Relation (H3) is a short form of the identity 



N 



N 



l\ sinh(A 7 — X a + 2rj) 



7=1 



N 



[] sinh(A a - Uj - 77) 



£ 

fc=i 



[] sinh(A 7 - v k + 77) 



7=1 



1 



ry . i / v Smh(A a - - T]) 



(45) 



2=1 



which may be proved, e.g., by induction, see also p5] . Relation ([43|) may be considered 
as the system of linear equations on the coefficients with the right hand side formed by 
the vector ^(Aq,), a = 1, . . . , N. Therefore, by Cramer's rule we obtain the connection 
formula for the determinants 

S-jv(Ai) 0(Ai,z/ 2 ) ... 0(Ai,i/jv) 
#jv(A 2 ) 0(A 2 ,z/ 2 ) ••• 0(A 2 ,^jv) 



det]yZ = — 



(46) 



JV 



9n{^n) 0(Atv, z/ 2 ) ... 0(Ajv,^jv) 

where det^i^ is given by formula (|2T]) . It is the transformed form of the determinant of 
the matrix Z which can be developed by the first column. Hence, one gets 

1 N 

det N Z = — ^(-lr^AS?.! (47) 
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V AT-1 



0(Ai,^ 2 ) 

0(A/9_i,l*j) 



0(Ai,Z/at) 

<f>(\p-h v n) 



(4* 



N-1 



(f)(X N , v N ) 

Substituting determinant representation (pi]) for the partition functions and Zjv-i 
in the left and right hand sides of relation Q5§D respectively, and cancelling the resulting 
common factors, one gets exactly equation fl47|). This proves that the determinant 
representation ( pip is the solution of recursion relation ( |5P| ) for the partition function. 

Let us turn now to the boundary correlation functions. Reduction formulae (p4T) 
and (|38|) , obtained in the previous Section, express them through the sum over partition 
functions of the models on (N — 1) x (N — 1) square sublattices. Using determinant 
representation fl2l"D for the partition function one can obtain then the determinant 
representations for the boundary correlation functions. 

Consider the correlation function H^. Substituting expression ([H]) for Z/v-i in 



reduction formula ( p4|) and extracting a general multiplier out of the sum over (3 one 
gets the following expression 



M-1 



N 



= Z N Y sinh 2rj J^J sinh(A Q — V\ 



rf) ] [ sinh(A c 

a=M+l 



vi + rf) 



X 



N N 

[] EI sinh(A Q - v k + rf) sinh(A Q - v k - rf) 

a=l k=2 

1] sinh(A/3 - A Q ) [I sinh(z/ fc - Uj) 

l<a<3<N 2<k<j<N 
M-1 N 

m 11 sinh(A 7 - A/3 + 2rf) J] sinh(A 7 - \p) 

7=1 



E 

/9=1 



iV 



(49) 



[] sinh(A /3 -u k -ri) 



k=2 



where the quantity ^^-1 ls defined by expression (f48|). Clearly, the sum over (3 in 



formula ( f4*9"D is nothing but the determinant of some N x N matrix developed by the 
first column. Only M first entries in this column are not equal to zero so there are M 
terms in the sum. Taking this into account and using expression ( |2"T| ) for Zn we finally 
obtain the following determinant representation for the correlation function Hj^: 



N 



H 



(Af) 
JV 



sinh2?7 Yl sinh(z/i 

k=2 



Vk) 



M N 

Y[ sinh(A a — ux + rf) ] [ sinh(A a — V\ 

a=l a=M 



v) 



detN^H 



where the entries of the N x N matrix Ji are given by 

'Hal = h M (Xa), H a k = 0(A a , Vk)-, k 



N. 



(50) 



(51) 
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The function h M (X) is equal to 

A/-1 N 

I] sinh(A 7 - A + 2r/) f] sinh(A 7 - A) 

M A ) = ^ . (52) 

n sinh(A -v k — rf) 

k=2 

The points Xm+i, ■ ■ ■ ,Xn are zeros of the function /ia/(A), hence, the last N — M entries 
in the first column of the matrix 7i are equal to zero. Note that the pref actor in (jjOp 
is, in fact, equal to 4>(Xm, ^i)/^a/(^i + ??)• 

Consider now the correlation function G^f\ Substituting expression ([H]) for Z N _i 
in reduction formula ([38]) we arrive at the expression for the correlation function 
as the sum of A^l r This expression is similar to (]4£f ) for the function Hjff\ Finally, 
we obtain the following determinant representation for the correlation function G^f . 

N 

II sinhOi -v k ) 

r (M) _ k=2 aetjyC/ ( zo\ 

N ~ M N det N Z' 1 ' 

Yl sinh(A Q - V\ + rf) Yl sinh(A Q — Vi—rj) 

a=l a =M+l 

where the entries of the N x N matrix Q are given by 

Gai = 9m{K), Gak = <f>{X a ,v k ), k = 2, . . . , N. (54) 
The function <?m(A) is equal to 

M N 

n sinh(A 7 - A + 2rj) Yl sinh(A 7 - A) 

/ \ \ 7=1 j=M+l , s 

^m(A) = - • (55) 

J] sinh(A -v k —rj) 

k=l 

At M = N this function is exactly the function ^jv(A) defined by formula fl42|) . The 
points Ajvf+ij • • • , Ajv are zeros of the function g M (A), hence, the last N — M entries in 
the first column of the matrix Q are equal to zero. By the direct check it is easy to get 
convinced that representations (p0|) and ( j53|) satisfy relations fl26|) and (p7|). 

Determinant representations (|50| ) and (^) for the boundary correlation functions 
and are the main results of the present paper. 



5. The free fermion case 



The "free fermion" condition implies the following restriction on the vertex weights: 

a 2 (A Q ,z/ fc ) + b 2 (A a ,z/ fc ) = c 2 {X a ,u k ), a, k = 1, . . . , N. (56) 

This equality is satisfied if we put r\ = if in ([I]). It is convenient to change the variables 
X a — > iX a , v k — > iz/fc, and after the rescaling a — > — ia, b — > — ib, c — > — ic one gets the 
following parametrization of the vertex weights: 

a(A a , v k ) = sin (A Q - v k + | ) , 

b(A OJ v k ) = sin (A Q - v k - f ) , (57) 
c(A a ,i/ fc ) = 1. 
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Since under the condition 77 = i| the determinant in equation becomes the Cauchy 
determinant, the partition function for the free fermion case can be evaluated explicitly: 

Zjv({A a }^ =1 ; {vk}k=\) = ] f cos(A Q - Xp) || cos(v k -Vj). (58) 

l<a</3<Af l<k<j<N 

For the correlation functions ffjy and G^f of the model in the free fermion case it 



is worth to use directly representations (0) and fl38|). Substituting formula (|58| ) in 
representation (|33|), one gets for the correlation function 

M-l N 

J! sin(A Q -vx-\) n sin(A a - v x + f ) 

■pj(M) _ a=l a=M+l 



N ~ N 

Y[ cos(z/j — V\ 
i=i 

iV 

M EI S in (V - V 3 + M 

X 



-s-i n oos(A a - A,) « Sm(A " A ' J) 



«=M 



The similar formula can be easily written down for the correlation function G^f\ 

Consider the homogeneous limit of the correlation functions in the free fermion 
case. The homogeneous model is obtained by putting the variables in each set {X a }^ =1 
and {vk]k=i t° be equal: 

Xi — . . . — A at = A, ui — . . . — u N = v. (60) 

Since now all the vertex weights depend only on the difference X — v, without loss of 
generality one can take v = —j. Thus, we have the homogeneous model with the 
following parametrization of the vertex weights 

a(A) = cosA, b(A) = sinA, c(A) = 1. (61) 

Note that the partition function of the homogeneous model is independent of A and 
is equal to one, Zn = 1, see (|58|). 

The substitution z/^ = v = — f in expression (1591) leads us to the following 



intermediate expression for when all are equal while all X a are still different 

M-l N M JV _ 1 . M 

<> = n - n - a. x e /° s X/> n -t^- m 

a=Af 



The problem now is to obtain the limit of expression (|62|) when all X a tend to the same 
value. For this purpose it is convenient to rewrite fl6"2|) in terms of rational functions 
instead of trigonometric ones 

M-l M N 1 M 1 

n-En^n^-. <«*> 



Hj* = + 



/3=1 a=M - ■ -»««/> t =1 M « - «/» 
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where 

u a = tan\ a , a — l,...,N. (64) 

Therefore, one should take the homogeneous limit in the set {u a }^ =1 : 

u a — ► u, a — l,...,N. (65) 

For the subset {u a }^ =M we may simply put um = %+i — . . . — Un — u, while for the 
subset {tiajfj"/ we parametrize u a as 

u a = u - (M - a)e, a = 1, . . . , M - 1, (66) 

and the homogeneous limit corresponds then to the case of vanishing e. Since the 
pref actor of the sum over j3 in ( |63| ) is regular as e — > 0, one may put e = in this 



perfactor. Thus, for one gets 



M - M 



h {m) = (i+ u 2 )u m - 1 v - fr 1 

N K ' ^[l + u 2 - u{M - p)e] N ~ M + 1 e M - 1 JJa-|3 



M-l , „ \Af— 1 



Taking into account that 



M-l 



TT 1 ("If (-1)' (M-l\ 

l\a-/3 /3\(M-f3-l)\ (M-l)\\ ) 1 
one obtains 



H 



(M)_ (-1)^(1 + ^ 2 ) ^ (-l)f /M-l 

(M - 1)! u N ~ 2M + 2 ^ [(1 + m 2 )^- 1 - e(3] N - M+1 a™- 1 \ (3 



The sum over (3 in (|69| ) in the limit e — > becomes exactly the (M — l)-th derivative 
of a function /(#) with respect to a variable z 

M-l 



^ (-iyy(z-/fe) /m-i\ d^- 1 

feL i*r=i ^ £ J = di^ /W ' (?0) 

/3=0 ' 7 



where 



/0) = ..v-m-m > 2 = ^— ■ ( 71 ) 



z iV-M+l : M 



Hence, for the sum over /5 in (|69|) in the limit e — > one gets 

r v til! V (-l)^(iV-l)! / M \" 

e ^ ^ [(I + m 2 )^- 1 -^-^ 1 ^- 1 V /5 / (N-M)\ \l + v?J ' l J 



Finally, substituting expression (|72|) in ( p9|) and taking into account that u = tan A 
we obtain the following expression for the correlation function of the homogeneous 
six- vertex model in the free fermion case (]61f) : 

rj(M) f N - l\ , 2 \ N—M ( . 2 \\M-1 f „n\ 

= \M — 1 J ^ C ° S 7 ^ ' ' ^ ' 
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Similarly, for the correlation function Gjf one gets 

M 

■I) _ ■ 

'N 



: e£;;)ka)»k))» (74) 

K=l \ ' 

Obviously, connection formulae fl26| ) and fl27|) are fulfilled by (|73|) and ([741). 

It is clear from the obtained results that the boundary correlation functions 
significantly depend on both A and M even in the simplest free fermion case. This 
dependence of the function G^ is of especial interest in the thermodynamic limit 
when both N and M go to infinity so that the variable x = M/N runs through the 
interval < x < 1. Let us denote the function G^f in this limit as G(x). Representing 
( [71]) in the form 

=(^~]) (sin 2 A) M - 1 (cos 2 X) N ~ M a F a (l, —M + 1; iV — M + 1; - cot 2 A) (75) 



and using the integral representation for the hypergeometric function 2-F1, one gets 
G ^ " (iV-M-l)!(M-l)! (8m A) (C ° S A) 



1 

X 



I dt{l-t) N ~ M - 1 (l + tcot 2 A) M ~\ M = l,...,iV-l. (76) 







Applying the steepest descent method to representation Q7E|) we obtain that G(x) is the 
Heaviside step function: 

r 1 if ^ > 

G(x) =6(x- sin 2 A) , 0(f) = \ \ if f = . (77) 

[0 if i < 

This result means that the arrows at the boundary column are ordered (frozen) in the 
thermodynamic limit: all arrows are pointing down above the point x = sin 2 A, while 
below all of them are pointing up. 



6. Determinant representations in the homogeneous limit 

In this Section the homogeneous limit of determinant representations ([50]) and (jj3|) 
for the boundary correlation functions is considered in the general case. To obtain 
the homogeneous model one should put all X a equal to A and all v k equal to v in the 
inhomogeneous model (|l|). Without loss of generatity one may put also v — 0. Hence, 
the vertex weights of the homogeneous six-vertex model are given by 

a(A) = sinh(A + r]), b(A) = sinh(A - 77), c(A) = sinh 2rj. (78) 

The procedure of taking the homogeneous limit for the partition function have 
been elaborated in the papers ]T7|, [I8| . In this limit, the singularities coming from the 
denominator of expression (|2l|) are cancelled by the zeroes coming from the determinant 
since then all rows and columns of the matrix Z tend to each other. 
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It can be shown by Taylor expansion of the entries of the matrix Z that the partition 
function of the homogeneous model is expressed through the double Wronskian 

det Z 

%am)a* = , » „ , Z _ 9 0(A), (79) 



^ 2 7V-1 ' V—nom/QK 



[0(A)] iv n« 



2 



n=l 



where 0(A) = 0(A, 0), namely, 

sinh 2?7 , . 

<t> A = — TT^ r— 77^ V- 80 

smn(A + 7]) sinn (A — f]) 

Since the correlation functions have been expressed in equations (150) and (EBI) through 



determinants, one may apply the approach described in detail in JT8[ to obtain the 
correlation functions of the homogeneous model. 

To apply the procedure given in [18| with minimal modifications, it is convenient 
to consider the function Gpf instead of G^f \ where 

G N (Ai, . . . , An; ui, . . . , vn) — G N (Xn, ■ ■ ■ , Ai; vn-, ■ ■ ■ > V\ ). (81) 

Clearly, these functions are equal in the homogeneous limit. Comparing equations ( j53|) 
and (|81|) one can see that the function G^ may be written as 

AT-1 

, . II sinh^iv -v k ) ~ 
f,M k=i det 7v Q 

N ~ n n-m det 7v Z ' 1 } 

[] sinh(A a - u N + rj) J] sinh(A Q -v N -rj) 

a=N-M+l a=l 

where 

Qak = 0(A«, Vk), GaN=9M(Xa), k = 1, . . . , N — 1, (83) 

and the function #m(A) is defined as 

N N-M 

n sinh(A 7 - A + 2rf) f] sinh(A 7 - A) 

- /w 1=N-M+1 7=1 , n .. 

9mW = N • (84) 

n sinh(A -u k -T]) 

k=l 

The matrix Q differs from the matrix Z by the last column, and the first N — M entries 
in the last column of the matrix Q are equal to zero, Qni — ■ ■ ■ = Qnn-m = 0. 

Now, the homogeneous limit in the set {A a }^ =1 can be easily found following 
[p"8| . Representing the differences A 7 — \ a in the expression for the quantity <?m(A q ) 
as (A 7 — A) — (A a — A), one may consider the differencies X a — A as independent variables. 
In the limit X a — > A, a = 1, . . . , N these variables tend to zero and it can be proved 
that the entries in the last column of the matrix Z, after successive subtractions of the 
rows become the coefficients in Taylor expansion of the function 

m = ( _ 1)y (sinh e )^sinh" (e -2,) ^ 
\ [ sinh(e + A — Uj — rj) 
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at the point e = 0. This solves the problem of taking the homogeneous limit in the set 
{A a }^ =1 . The homogeneous limit in the set {vk}k=i can be taken in the same manner as 
for the partition function, with the only difference that there are no subtractions from 
the last column; one should simply put all equal to v = in the entries of the last 
column. 

As a result, one obtains the following determinant representation for the correlation 
function G^f in the homogeneous limit: 

r (M) _ (N - 1)! detjyg hom 

N " [sinh(A + T7)] M [sinh(A - rj)] N - M det N Z hom ' { ' 



where 



{Qhora)ak (-2^hom)afc; k 1> • • • j -A 1, 

_ d"" 1 J (sinh5) Ar - M [sinh(e - 2t])} m \ 



(87) 

£ = 

For the correlation function ifjy in the homogeneous limit the following expression is 
valid 

H (M) (N - l)!sinh2r] det N H hom 

N = [sinh(A + 77)] M [sinh(A - v )] N - M + l det N Z hom ' 1 j 



where 



(^hom)ofc (-2^hom)«fc; ^ !)•••> N 1, 

d"' 1 / (sinhg) jV - M [sinh( £ - 2 V )] M - 1 \ 
(H hom ) aN - | [sinh(e + A - V )] N -i J 



(89) 



£=0 



Formulae (pq ) and (|88|) generalize representation (|79| ) for the partition function of the 
homogeneous model, and they may be used for investigation of the correlation functions 
in the thermodynamic limit. 

7. Conclusion 

In the present paper the determinant representation for the one-point boundary 
correlation functions of the six- vertex model with the domain wall boundary conditions 
is derived. For the correlation functions Hff fl25f) and G^ (|24]) we have obtained 



representations fl50|) and (|53|) , which generalize the determinant representation for the 
partition function Z^. In these formulae the correlation functions are expressed through 
the determinants of the N x N matrices 7i and Q respectively. The matrices 7i and 



Q differ from the matrix Z appearing in representation fl2T|) for the partition function 
by the first column only. The derivation of representations ([50]) and (^) is based on 
reduction formulae (p4|) and (|38| ) which have been obtained in Section 3 exclusively by 
means of the algebra of operators entering the monodromy matrix. 

In the free fermion case studied in Section 5 we have obtained explicit (determinant- 



free) formulae ([F^) and ( |74D for the boundary correlation functions and G N of 

(A/' 
N 



the homogeneous model. In the thermodynamic limit N, M — > oo the function G^ M 
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describing the spontaneous polarization at the boundary turns into the function G(x) 
with x = M/N, < x < 1. We have found that G(x) is just the Heaviside step function 
(\n\). The emergence of the step function for the spontaneous polarization indicates the 
"freezing" of the arrows at the boundary in the thermodynamic limit. On the other 
hand, from the results of the paper ]2(J we have found numerically that at the ice point, 
a = b = c = 1, the function G(x) exhibits the similar behaviour, G(x) = 8(x — 1/2). 
The obtained determinant formulae (|86|) and (|88l) may be used for investigation of 
the boundary correlations in the thermodynamic limit in the homogeneous model with 
arbitrary values of the vertex weights. 

The step function behaviour of the boundary spontaneous polarization indicate the 



existence of the analogue of the arctic circle theorem [21], p2j| . To be more precise, we 



expect that for the values — 1 < cosh 27/ < 1 the arrows are "frozen" at the corners of 
the grid, while inside the grid there is a region of "disorder" . To obtain the shape of the 
"disordered" region of the grid it is necessary to obtain an appropriate expression for 
the spontaneous polarization not only at the boundary but also at the arbitrary point of 
the lattice. Hence, the problem of calculation of the correlation functions of the model 
deserves further investigation. We hope that the approach described in Sections 3 and 
4 may appear to be productive in the derivation of the correlation functions outside of 
the boundary. 
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